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;> . Abstract 

l/^ , We are studying the effect of converting caught targets into new chasers in the context of the 

o 

recently proposed 'group chase and escape' problem. Numerical simulations have shown that this 
f~^ ' conversion can substantially reduce the lifetimes of the targets when a large number of them are 

initially present. At the same time, it also leads to a non-monotonic dependence on the initial 
number of targets, resulting in the existence of a maximum lifetime. As a counter effect for this 
5_^ ■ conversion, we further introduce self- multiplying abilities to the targets. We found that the longest 

lifetime exists when suitable combination of these two effects is created. 
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I. INTRODUCTION 

The collective motions of interacting entities have recently attracted a lot of attention. 
They include groups such as molecular motors, ants, fish, birds, pedestrians, and vehicles. 
They are generally termed self-driven particles (SDPs). The universal phenomena in these 
systems, such as the phase transitions and metastable states, have been discovered and 
investigated [1-3]. 

'Chase and Escape', often observed in nature, has long been studied in mathematical 
literature. Even the simplest of cases, a two particle system composed of one chaser and 
one target, posed the challenging mathematical problem of analytically describing their 
trajectories [4-6]. This chase and escape problem can be apphed to SDPs and also systems 
that show chemotaxis [7]. Systems with a number of chasers and one target have been 
modeled and analyzed [8,9]. Systems with several chasers and targets were investigated in 
the fields of game theory, mechanical engineering, and multi-agent problems [10,11]. This 
topic has recently been further extended to take into consideration cases where large numbers 
of chasers and targets exist, which is called "group chase and escape" [12]. In this model, 
the chasers pursue the nearest targets, while the targets try to escape from be caught 
by the nearest chasers to avoid removal. We should note that in a study of collective 
motion, pursuit-and-escape interactions were also investigated by using Brownian particles 
not grouped separately as chasers and escapees [13]. In contrast, in group chase and escape, 
particles are grouped separately. Some of the properties associated with such group chase 
and escape were found. For example, spontaneous self-organizations are formed for both the 
chasers and targets without any internal communications within each group. In addition, 
given the initial number of targets, the optimal number of chasers can be identified by 
introducing the running cost. We have extended this model of group chase and escape in 
this paper by posing to ourselves the following question, "How do the dynamics change if 
the caught targets are converted to chasers, instead of removing them?" This conversion can 
be related to the spread of rabies [14] by interpreting that the chasers are infected animals 
and the targets are susceptible ones. The targets caught by the chasers become infected, 
i.e., they are converted to chasers. To answer the above question, we investigated how the 
lifetimes of targets change as we vary the rate of conversion. This article is composed as 
follows. Our model with the given conversions is introduced in section [TTl Then, we discuss 



the numerical simulations and their results in section UTTl Section HVl contains our conclusive 
discussions. 

II. MODEL 

We took into consideration two kinds of particles: chasers and targets. N^ chasers and 
Nj, targets are initially randomly placed on a L^. x Ly discrete square lattice with periodic 
boundaries. We define the discrete time update for the chasers and targets. At each time 
step, the chasers and targets move according to the following rules. Each chaser at {xc, He) 
tries to catch its nearest target. The nearest target is a target placed at a minimum distance 
from the chaser, where the distance between two points (xi,?/i) and (x2,|/2) is defined by 
using the following equation, 

^{min[a;i - X2, L^ - (xi - X2)]}^ + {min[yi - ^2, Ly - {yi - 1/2)]}^ 

If there are two or more targets with the same minimum distance from a chaser, one of them 
is randomly chosen as its nearest target. The position of the chosen target is denoted by 

{XCn,yCn)- 

The behavior of each chaser changes whether the chosen target is located next to it or 
not (i.e., the distance is one or larger). If the distance is larger than one, it hops one of 
its four neighboring cells to decrease the distance, as shown in Figure [H If there is more 
than one such cell, one of them is randomly chosen using equal probabilities. However, we 
also assume exclusive volume effects. So, the chaser can hop only if the cell is not occupied 
by another chaser. If the cell is occupied, it stays at its current position. On the other 
hand, if the target is located next to the chaser, it is caught. In the previous study [12], 
this catch event resulted in the removal of the target. In this paper, we introduce the 
new additional possibility of converting the target into a new chaser, as shown in Fig. [21 
'Remove' or 'Convert' actions take place with the probabilities 1 — Py and Py, respectively. 
If a remove action is chosen, the chaser hops to the position of its nearest target and the 
target disappears. If a convert action is chosen, the chaser does not move and the nearest 
target is converted into a new chaser. We note that Py = corresponds to a situation where 
all the caught targets are removed leading to the original model [12]. Py = 1 represents 
the case where all the caught targets become chasers. The movements of the targets are 
as follows. A target at {xtiVt) tries to hop one cell to escape from its nearest chaser at 



{xTn,yTn), and if there are multiple chasers with the minimum distance to the target, one 
chaser is randomly chosen as the nearest chaser. Similar to the chasers' rule, the target hops 
one of its four neighboring cells to increase the distance, as shown in Fig. [T] If there is 
more than one such cell, one of them is randomly chosen using equal probabilities. We also 
assume exclusive volume effects for the targets. We use the following overall discrete-time 
updating schedule. In each time step, all the chasers are moved in a random order first, and 
then, all the targets are updated using a random order. 

III. SIMULATIONS 

We investigated the entire lifetime T and typical lifetime r for various conversion prob- 
abilities, Pv, using numerical simulations to evaluate the effect of the conversion rule. T is 
defined as the period from the start time t = to the time when the last target is caught. 
T is defined as E^i i{Nt^ - N^)/N^. 

We plotted T and r versus A^"^ with Py G {0,0.01,0.02,0.05,0.1,0.2,0.5,1}, in Figs. 
[3] and HI respectively. The following conditions were given. The size of the system is 
Lx = Ly = 100, the initial number of chasers is A^^ = 100, and the results are averaged over 
1000 samples. In addition to the chasing and escaping rule defined in section [TTl we also 
present the results for a case in which all the chasers and targets obey random walks for a 
comparison. 

Maxima T exist for the plots in the Py > 0.01 range for both the chasing and escaping 
rule and the random walking rule. However, there is no observed maximum T for Py = 0. As 
for r, the maxima r exists for all the Py for the chasing and escaping rule. For the random 
walking rule, however, r stays almost constant at Py = and monotonically decreases with 
A^° for Pv > 0.01. 

The existence of a maximum T can be explained as follows. In the original group chase 
and escape model {Py = 0) [12], T increases as A^^^. If A^^ is small, the conversion Py > 
would not be relevant. However, as A'^j. increases, the chasers frequently catch targets and 
convert them to new chasers, which results in the fast spread of a chain of conversions. This 
spread would drastically decrease the entire lifetime T. Thus, these two factors lead to the 
maximum T as a function of A^^^. This would also be applied to the existence of a maximum 
T for the chasing and escaping rule. The latter factor also causes the monotonic decrease 



r for the random walking rule. Note that a maximum r in Py = was reported in [12]. 
Furthermore, it is noteworthy that T with Nj, = 1 is larger than T with A^j. = 5000 for a 
large Py for the random walking case. 

We also plot T and r as a function of A^^. The initial number of targets is fixed at 
A''^ = 100, the degree of conversion is given as Py G {0, 0.2, 0.5, 1} and the other conditions 
are the same in Figs. [3] and HI 

The T and r results versus A*"^ are shown in Figs. [5] and [6l respectively. For the chasing 
and escaping cases, a T with various Py merges around N^ = 40 and a r also merges around 
Nq = 40. The lines do not depend on Py above the critical points. Meanwhile, there is 
a spread in the slope below the points. Still, in the case of chase and escape, we observed 
clear changes in the slopes around the points, while in the random walk case they smoothly 
change. 

Such clear changes for the chasing and escaping rule were reported for Py = [12]. This 
difference reflects the sudden changes in effectiveness of catching the targets. That means 
there is a critical number of chasers, beyond which adding more chasers is not as effective. 
Our results indicate that this critical number also exists even with a conversion, Py > 0. 
In addition, this understanding is consistent with the result showing that the critical point 
does not depend on Py, but is dictated by a large number of chasers. 

In the above simulations, the targets monotonically decrease due to the catch-up and 
are bound to extinction. In addition, the conversions from targets to chasers are more 
advantageous for the chasers. Here, we consider giving the targets the ability to proliferate 
in order to resist against extinction. The details for this new rule are given as follows. When 
each target hops to one of its neighboring sites in a single time step, a new target arises 
on the original site with the probability Pt- Pt is a parameter for the natural increase. 
Pt = is the case in which the targets do not increase, while Py = 1 is the case in which the 
targets always increase when they move. Note that this increasing rule takes into account 
the exclusive volume effect. When the site is filled with a target, a new target does not arise 
due to the shortage of space. Note also that the new target stays in the time step. 

We plot the relationships between the lifetime T and Pt for various Py. The conditions 
for this are given as follows. The size of the system is L^ = Ly = 50, the initial numbers are 
A^^ = 10 and A^j. = 50, and the simulations are performed for 1000 samples. The results are 
shown in Figs. [7] and [S] for the chasing and escaping and random walking rules, respectively. 



In both figures, there is a maximum of T for a small Py- 

The mechanism for producing such a maximum T can be explained as follows. If Pt 
is very small, all the targets are caught by the chasers before they can proliferate. As Pt 
increases, the targets have more chance to increase, which leads to an increase in T. With 
Prp increasing, however, the converted chasers make T smaller. In this case, the targets 
initially explosively increase, and afterward a number of targets are caught and converted 
to new chasers. Due to this increase in converted chasers, the catching rates increase, and 
finally the number of targets start to very quickly decrease. The existence of a maximum 
T for a small Py suggests that the proliferation rate of the targets should not be too high 
for longer survival. On the other hand, if Py is large, the conversion ability of the chasers 
overrides the targets' proliferation abilities, which obscures the maximum. 

IV. CONCLUSIVE DISCUSSIONS 

We proposed group chasing and escaping models with a conversion rule from targets 
to chasers. We found, through numerical simulations, that the conversion rule provide the 
optimum initial number of targets to produce the maximum lifetime T for the targets. While 
this maximum is also observed in the random walk rule, the conversion does not affect the 
qualitative difference with respect to the changes in the initial number of chasers. We also 
performed simulations using the effect of self-proliferation of the targets. The simulation 
results again exhibited the maximum lifetime at the appropriate level of proliferation when 
the conversion rate was moderate. 

We will take the followings into consideration as future works. First, we need to explore 
whether real data exhibit the qualitative behavior of the model. Second, we want to take into 
consideration a time delay element, which in reality may be relevant in infective processes, 
such as a conversion delay. Third, we would also like to extend our model to include a 
variety of target resistances, which in turn would be beneficial for developing more efficient 
strategies for the chasers. 
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FIG. 1: Directions of hopping for chasers (red) and targets (green). 
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FIG. 2: A caught target (green) is converted to a new chaser (red) with probabihty Py- 
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FIG. 3: Relationships between T and A'j. with various Py for the two cases, the chasing and 
escaping case (C&E) and the case of random walks of chasers and targets (RW). 



10 



200 
150 

100 
70 
50 

go 

20 

15 

10 









C&.E:P\/=0 

Ce.E:a'=0.01 

Ce.E:P|;=0.02 

Ce.E:P|;=0.05 

Ce.E:P|;=0.1 

Ce.E:P|;=0.2 

Ce.E:P|;=0.5 

Ce.E:P|/=1 



RW:PV=0 

RW:PV=0.01 

RW:PV=0.02 

RW:PV=0.05 

RW:PV=0.1 

RW:PV=0.2 

RW:PV=0.5 

RW:PV=1 




10 



100 




1000 



40^ 



A/: 



7 



FIG. 4: Relationships between r and A'^j. with various Py 
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FIG. 5: Relationships between T and N(^ with various Py- The notations C&E and RW correspond 
to the chasing and escaping case and the case of random walks of chasers and targets, respectively. 
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FIG. 6: Relationships between r and N9^ with various P] 
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FIG. 7: Relationships between T and Pt with various Py in the chasing and escaping case. 
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FIG. 8: Relationships between T and Pt with various Py in the case of random walks of chasers 
and targets. 
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